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ABSTRACT 

We present an analysis of a series of four consecutive Chandra high resolution transmission grat- 
ings observations, amounting to a total of 150 ks, of the Be X-ray source HDl 19682 {— IWGA 
J1346.5— 6255), a member of the new class of 7 Gas analogs. The Chandra lightcurve shows signifi- 
cant brightness variations on timescales of hours. However, the spectral distribution appears rather 
stable within each observation and during the whole campaign. A detailed analysis is not able to 
detect any coherent pulsation up to a frequency of 0.05 Hz. The Chandra HETG spectrum seems 
to be devoid of any strong emission line, including Fe Ka fluorescence. The continuum is well de- 
scribed with the addition of two collisionally ionized plasmas of temperatures kT w 15 keV and 0.2 
keV respectively, by the apec model. Models using photoionized plasma components (mekal) or non 
thermal components (powerlaw) give poorer fits, giving support to the pure thermal scenario. These 
two components are absorbed by a single column with N}i = {0.20~^q q^) x 10^^ cm~^ compatible with 
the interstellar value. We conclude that HD119682 can be regarded as a pole-on 7 Gas analog. 
Subject headings: stars: individual (HD119682) - X-rays: binaries 



1. INTRODUCTION 

HD119682 (= IWGA J1346.5-6255) is a Be star which 
presents copious X-ray emission. A member of the open 
cluster NGG 5281, located at a distance of ^ 1.3 kpc, 
this Be star displays an X-ray luminosity in the 0.5-10 
keV energy band of the order of Lx ^ 10"^^ ergs~^, low 
when compared with the typical luminosities of Be X- 
ray binaries (Lx ~ iq^g-st ergs~^), two orders of mag- 
nitude lower than persistent Be X-ray binaries (Reig & 
Roche, 1999) and at least one order of magnitude higher 
than the expected X-ray emission from isolated OB stars 
(Berghofer et al. 1997). With a mass of M = 18 ± IMq 
and an age of 4 ± 1 Myr, HD119682 seems to be a blue 
straggler (Marco et al. 2009), the second (out of the 7 
known) gamma Gas analogs that shares this property. 

Rakowski et al. (2006; R06), using the superb spa- 
tial resolution of Chandra as well as Ha images from the 
Glay 6.5 m Telescope in the Magellan Observatory at Las 
Gampanas, unambiguously identified the X-ray source 
IWGA J1346. 5-6255 with the OB star HD119682. Us- 
ing XMM-Newton PN+MOS data, these authors showed 
that the X-ray spectrum is well described by the addi- 
tion of two optically thin (raymond) thermal plasmas 
with kTi = 1.07 ± 0.16 keV and fcTa = 10.41^-3 keV 
respectively. Based on the similarities shared between 
HDl 19682 and 7 Gas (early type optical counterpart, 
thermal nature of the X-ray emission, log Lx /iboi = 
—5.4, etc) R06 classify the source as a 7 Gas analog, 
and argue against other X-ray production mechanisms 
like colliding winds (with an unseen companion) or emis- 
sion from a magnetically confined wind (given the very 
high temperature of the hot plasma) . 

Safi-Harb et al. (2007; S07), in a multifrequecy 
analysis, confirm the identification between IWGA 
J1346.5— 6255 and HD119682 and characterize the op- 
tical counterpart (B0.5Ve), definitely unrelating this 



source from the background Supernova Remnant G309.2- 
00.6. These authors find that the XMM-Newton 
PN-hMOS data can be satisfactorily described by the 
addition of two optically thin thermal plasmas, using 
the mekal model, with kTi = 1.7 ± 0.3 keV and kT2 = 
13.0^2 4 keV respectively, while the fit to the Chandra 
spectrum, extracted from the S3 AGIS GGD chip, deliv- 
ers temperatures of kTi = 0.95^1^^4 keV and fcTa > 14 
keV respectively. Alternatively, the hotter thermal com- 
ponent can be described by a powerlaw. This degeneracy 
has also been found in 7 Gas (Smith et al. 2004, Lopes de 
Oliveira et al. 2011) and HD110432 (Lopes de Oliveira 
2007, Torrejon et al. 2012). S07 present also hints of 
a possible pulsation with P ~ 1500 s which could be 
caused by the spin of a NS or a WD. 

The origin of the X-ray emission in 7 Gas analogs is 
still a mystery. Two hypotheses have been put forward: 

• The X-ray binary scenario in which a White Dwarf 
(WD) accretes matter from the Be donor (Haberl, 
1995; Kubo et al., 1998; Owens et al., 1999 for the 
case of 7 Gas; Torrejon & Orr, 2001 for HD110432). 
The low potential well of the WD, compared with 
that of a NS, accounts naturally for the lower lumi- 
nosity and the thermal nature of the X-ray emis- 
sion. This explanation faces several problems. On 
one hand, the X-ray emission shows no sign of any 
orbital modulation or spin pulse. On the other 
hand, the progenitor of the putative WD should 
have been still more massive than the early B com- 
panion (M ^ 2OM0). From the evolutionary per- 
spective, it is challenging to produce a WD from 
such a massive star; naturally, they would explode 
as a supernova leaving behind a NS, as is the case 
for the Be X-ray binaries. In the case of HD119682 
no SN explosion seems to have occurred (Marco et 
al. 2009) which argues against the NS companion. 
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TABLE 1 

Journal of Chandra HETG observations. 



ObsID 


Date 
(dd-mm-yyyy) 


(ks) 


Rate" 
(e/s) 


8929 


17-12-2008 00:20:10 


28.9 


0.039 


10835 


19-12-2008 09:32:42 


29.7 


0.031 


10834 


20-12-2008 15:45:04 


59.3 


0.032 


10836 


21-12-2008 15:47:03 


29.9 


0.035 



^ Zcroth order 



• The standalone star scenario in which the X-rays 
are produced by the interaction of the star's mag- 
netic field with the circumstellar disk (Robinson 
et al. 2002, Smith & Robinson 2003, Smith et 
al. 2004). In this scenario, magneto rotational 
instabilities in the inner part of the circumstel- 
lar disk around the Be star (close to the Keple- 
rian co-rotation radius) enhance a seed magnetic 
field through a disk dynamo mechanism (Balbus & 
Hawley, 1998). The differential rotation of the disk, 
stretch and sever the entangled field lines acceler- 
ating particles which, upon impact on the star's 
surface, heat the plasma up to MK temperatures. 
The finding that the X-ray emission seems to be 
produced close to the surface of the Be star sup- 
ports this scenario (Smith et al. 2012, Torrejon et 
al. 2012). The details of this mechanism, though, 
are far from clear. Why other early type Be stars 
with large circumstellar disks do not produce X-ray 
emission at all is not explained. 

In any case, to understand the true nature of these 
systems will have important consequences in our picture 
of the structure and evolution of early type stars and/or 
how binary systems evolve. 

In this work we present an analysis of a series of four 
consecutive Chandra High Energy Transmission Gratings 
(HETG) observations of the Be star HD119682 allow- 
ing us to study the spectrum at the highest resolution 
available today. The four observations amount to a total 
of ^ 150 ks, spread over five days, allowing a sensitive 
search for pulsations in a wide range of frequencies. 

2. OBSERVATIONS 

We have processed and analyzed the Chandra data 
from the ObsIDs 8929, 10835, 10834 and 10836 avail- 
able at the public archive (PI C.E. Rakowski). In Ta- 
ble [1] we show the journal of observations. The Chan- 
dra High Energy Transmission Gratings spectrometer 
(HETG; Canizares et al. 2005) observed the source dur- 
ing 150 ks, splitted in four intervals. Both HETG in- 
struments, the High Energy Gratings (HEG; 0.7-8 keV) 
and Medium Energy Gratings (MEG; 0.4-8 keV) were 
used for the analysis. After filtering the data, the 0.5 
to 7.5 keV spectral range showed enough signal-to-noise 
ratio for the scientific analysis. Both the spectra as well 
as the response files (arf and rmf) were extracted us- 
ing the CIAO software (v 4.4). First dispersion orders 
(to = ±1) were fitted simultaneously. The extracted 
spectra were binned to match the resolution of HEG 
(0.012 A FWHM) and MEG (0.023 A FWHM) and have 
a minimum S/N ratio of 7. The peak count rate in a sin- 
gle (MEG) gratings arm spectrum is w 0.008 cts/sec/A 



and, consequently, the spectra is not affected by pileupQ- 
For the lightcurve extractions we add data from all dis- 
persed orders (to = — 3, — 2, — 1, 1, 2, 3) from HEG and 
MEG. High dispersion orders (to = ±3, ±2) constitute a 
25 % of the total count rate. 

The analysis was performed with the Interac- 
tive Spectral Interpretati on System (ISIS) v 1.6.1-24 
(jHouck fc Denicolal |2000[) . The photoelectric absorp- 
tion has been treated with the tbnew absorption model 
(Wilms et al. 2012), which includes the most up to date 
photoelectric cross sections. 

3. TIMING ANALYSIS 

The four observations analyzed in this work span five 
different days with uninterrupted coverages ranging from 
30 ks to 60 ks. In total they represent 150 ks worth of 
data. This allows us to search for coherent pulsations 
over a wide range of frequencies with unprecedented de- 
tail. In the left panel of Fig. [T] we present the Chandra 
lightcurve binned to 2000 s bins. As can be seen, the 
source brightness is variable on timescales shorter than 
the duration of a single observation (ObsID). During the 
whole campaign, however, the source remains quite sta- 
ble. In order to search for coherent pulsations we have 
analyzed the lightcurve with a maximum time resolution 
of 10 s. Subsequently we have created lightcurves with 
larger time bins of 100, 250 and 500 s to reduce the error 
bars. In Fig. [l] right panel, we show, as an example, 
the normalized Fast Fourier Transform periodogram of 
the lightcurve using 250 s bins, using the Lomb - Scargle 
technique, for each ObsID. The upper plot corresponds 
to all four intervals together taking into account the du- 
ration of the time gaps between each observation. Safi- 
Harb et al. (2007) suggested the possibility of a ~ 1500 s 
pulsation in the lightcurve of HD119682. In our last Ob- 
sID (10386) there is indeed a peak at P = 1443±17s (and 
also other at 2350±50 s). The significance of this peak, as 
given by the False Alarm Probability (FAP), is low. FAP 
is defined as the probability that no periodic component 
is present in the data with such a period within the cur- 
rent frequency search parameters. In order to compute 
this significance we have run 10^ Monte Carlo simula- 
tions of the randomized Y-axis values of the lightcurve 
and performed the corresponding periodograms. Then 
we search for peaks higher than that in the original 
periodogram at any frequency. We have found higher 
peaks in 20 ± 3 % (two sigma standard error) of the pe- 
riodograms. Therefore, the False Alarm Probability is 
rather high, ^ 20% and, consequently, the significance of 
the peak low. Similar results are obtained irrespectively 
of the size of the time bin. However neither of these peaks 
are present in the previous data sets. Although some low 
significance peaks can be seen in each periodogram, none 
of them remain stable. In consequence, in the whole data 
set we do not find evidence of any significant period up 
to a frequency of 0.05 mH^. 

In order to study the spectral variations of the X- 
ray source we have extracted the Chandra lightcurves 
in the wavelength intervals 3-6 keV (H), 2-3 keV (M) 
and 0.5-2 keV (S), and have created the correspond- 

1 See The Chandra ABC Cuide to Pileup, v.2.2, 
http : / /cxc .harvard. edu/ciao/download/doc/pileup-abc .pdf 

The Nyquist frequeney eorresponding to the 10 s time bin. 
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Fig. 1. — Left: Chandra lightcurve, in the 2 - 25 A wavelength range, binned in 2000 s intervals. All lightcurves have an arbitrary time 
offset and, therefore, the time gaps are not to scale. Right; Power spectrum corresponding to the 250 s binned Chandra HETG lightcurve. 
The 1500 s pulsation {v ~ 6.7 X 10"'* Hz) is only detected in ObsID 10836 but is not detected neither in the rest of the ObsIDs nor in the 
combined ligthcurve. 

ing hardness ratios versus source intensity plots. In Fig. 
[21 we show the variation of the soft/high hardness ra- 
tio against the source intensity (0.5-6 keV) using 7000 s 
time bins. As can be seen, the source does not show any 
spectral correlation with brightness. The average val- 
ues of the hardness ratios are 0.43 ± 0.07 (ObsID 8929), 
0.48±0.07 (ObsID 10835), 0.53±0.06 (ObsID 10834) and 
0.50 ±0.06 (ObsID 10836). They are, therefore, compat- 
ible within errors. The differences between single values 
of the hardness ratio are not significant either. Individ- 
ual errors range between 0.15 and 0.20. In conclusion, 
while changes in brightness from bin to bin are, overall, 
significant the hardness ratios are consistent with a con- 
stant value. The source appears rather monochromatic 
on timescales from hours to days. 
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4. SPECTRAL ANALYSIS 

According to the previous analysis, the spectra of 
HD 11 9682 is fairly stable in a time scale of days de- 
spite the variations in the brightness. Therefore we have 
combined the four observations into a single spectrum to 
increase the signal to noise ratio (S/N). The resulting 
spectrum is shown in Fig. [3l 

In order to describe the continuum, models consist- 
ing of one component do not give satisfactory results. 
For example, an absorbed powerlaw gives a reduced chi 
square Xr = l-'^'l for 88 degrees of freedom (dof) while 
a single mekal gives xl = 1-64 (88 dof). The spectrum 
shows always a pronounced excess below 1 keV which 
can only be modeled when a second component to the 
model is added. The best fit is attained with the ad- 
dition of two optically thin collisionally ionized plasmas 
(apec). The parameters are shown in the fourth col- 
umn of Table [2] The low energy range is fitted by a 
plasma with kT = 0.20lEJ:j^^ keV and solar abundance. 




0.02 



0.03 0.04 
Counts/s 



0.05 



Fig. 2. — Chandra hardness-intensity plot using 7000 s time bins. 
Different colors represent the different ObsIDs, with the same color 
coding as in Fig. [l] While changes in brightness from bin to bin 
are, in general, significant the hardness ratios are consistent with 
a constant value. 

The high sensitivity of the Chandra gratings allow us to 
constrain this temperature with unprecedented accuracy. 
Although there are few energy bins below 1 keV, the as- 
sociated errors (blue crosses) are very small. The ob- 



tained plasma temperature kT2 = 0.20_ 



0.02 



keV is fully 



compatible with the temperature of the cool plasma in 
IID110432 derived from high resolution spectra (Torrejon 
et al. 2012) and significantly lower than that obtained 
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TABLE 2 

Model parameters for Chandra HETG data. 



Component 


Parameter 




Value 






PL+APEC 


2 APEC 






0. 15 Q 




powGrlciw 


Norm 

r 


00026+^-''??^ 
u.uuu<iu_o 0004 

'^■•^"-0.01 




apec 1 


Norm 
kT (keV) 




0.0018lHgo6 


apec 2 


Norm 
kT (keV) 


^ + 0003 


r, r./-vr.i-4-0 0001 

0-0005Io:ooo2 




Flux' 


1.78 


1.81 




X? (dof) 


1.48 (86) 


1.35 (86) 


^ In units of X 10^^ cm"^ 

^ Unabsorbcd 0.5-8 kcV fiux in 


units of X 10^^^ org 


s^^ cm^^ 
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Fig. 3.— The ChandraHETG unfolded spectrum of HD119682 
showing the relative contributions of the low and high temperature 
plasmas respectively. Also shown (in red) is the best fit model. 

by R06 {kT = 1.07 ± 0.16 keV; Raymond-Smith) and 
S07 (kT = 1.7 ± 0.3 keV; mekal). The absorption col- 
umn, however, remains essentially compatible (see be- 
low). Therefore, the differences in kT could be due to 
intrinsic long term variability of the source. This vari- 
ability in the cool plasma temperature seems also to be 
present in 7 Cas (Lopes de Oliveira et al. 2010) while 
the hot plasma temperature appears to be more stable. 

The high energy range requires a second, much hot- 
ter apec model with kT = 15.7^g'4 keV. This tempera- 
ture is compatible within the errors with that derived by 
R06 and S07. Alternatively, in agreement with previous 
works on 7 Cas candidates, the high energy range can 
be fit with a powerlaw with a photon index F = 1.36. 



The parameters are shown in the third column of Ta- 
ble [2] However, in terms of the two apec model is 
significantly better. Replacing the coUisionally ionized 
plasmas (apec) by photoionized plasmas (mekal), gives 
slightly worse fits. For example, the two mekal model 
gives Xr = 1-38 (86 dof) while powerlaw + mekal gives 
Xr = 1-49 (86 dof). In conclusion, the fully thermal 
model composed by two apec components is favored by 
the present data set. 

The most conspicuous characteristic of the spectrum 
is the lack of any discernible emission line. In particu- 
lar, no Fe Ka line is seen and even the fiux at 6 keV is 
below the 7 a level. There is flux up to 7 keV at the 3tT 
level; however, there is certainly not an Fe emission line 
even at 1 a. Finally, the residuals present some structure 
between 1 and 2 keV that could be due to the presence 
of emission lines. Indeed, this is the line rich region for 
highly ionized plasmas. The HETG spectra of bright 7 
Cas analogs display several strong emission lines in this 
region (Smith et al. 2004, Fig. 2; , for 7 Cas; Torrejon et 
al. 2012, Fig. 7 and 8, for HD110432). Strong lines of Ne 
X Lya (1.02 keV) and Si xiv Lya (2.00 keV) as weU as 
weaker lines of Fe xxiv (1.17 keV) and Mg xii Lya (1.47 
keV) might be expected in this area. However, none of 
these lines show up in the spectrum of HD119682 with 
EW upper limits of the order of 1 eV. 

The high sensitivity of Chandra HETG at low ener- 
gies also allows us to study the absorption column with 
unprecedented detail. The analysis of brighter mem- 
bers of this class have shown that the different spectral 
components are affected by different absorption columns 
(Smith et al. 2004, for 7 Cas; LO07 and Torrejon 
et al. 2012, for HD110432). Here we have tried the 
same strategy. The values of the resulting columns turn 
out to be identical, within the errors, for all models. 
Therefore a single photoelectric absorption satisfacto- 
rily describes the spectrum of HD119682. The measured 

with R06 and S07 



Q3) X 10 cm ^ is very low, in agreement 



5. DISCUSSION 

The spectrum of HD119682 differs from other 7 Cas 
analogs observed at high resolution mainly in that it 
shows no discernible emission lines. Particularly strik- 
ing is the fact that no Fe complex is observed. Despite 
this, the continuum spectral parameters of HD119682 
are very similar to those of 7 Cas (Smith et al. 2004, 
Fig. 2) and HD110432 (Torrejon et al. 2012). Consis- 
tently with these previous works, the correct description 
of the spectrum is attained with the addition of several 
thermal components. In the case of HD119682 only two 
are needed with characteristic temperatures in the ~ 2 
MK and ^ 150 MK range. Both components are mod- 
ified by a single absorption column ~ 0.20 x 10^^ 
cm~^. This column would imply an interstellar ex- 
cess E{B -V) ^ N^/6M X 10^1 = 0.29 j:!^:"^^ (Ryter 
199(13) in good agreement with the interstellar medium 
(ISM) absorption towards NGC 5281 (Moffat & Vogt 
1973, E{B - V) ^ 0.26 ± 0.03; Marco et al. 2009, 
E{B -V)^ 0.28 ±0.02). 

^ The errors represent the propagated uncertainty in the mea- 
sured 



Chandra observations of HD119682 



5 



HD119682 has been observed several times both in X- 
rays {XMM-Newton August 2001 and Chandra Decem- 
ber 2004, R06, SOT; December 2008, this work) and in 
the Ha Ught (June 2003 and February 2004; R06, S07) 
showing always a strong emission. Although the time 
coverage is still poor, the source appears to be persistent 
in both energy bands. Assuming that the circumstellar 
envelope (and hence the Ha emission) was also present at 
the epoch of our Chandra observation, the X-ray source 
was not obscured by the disk material. The same would 
be true for the observations analyzed in previous works 
(R06, SOT). The circumstellar disk (confined to the equa- 
torial plane) does not intercept the X-ray emission. In 
the framework of the accreting binary system, the com- 
pact object (the WD) should not be deeply embedded 
into the circumstellar disk; yet it should be accreting 
enough material to power persistently the X-ray source 
up to a non negligible Lx ~ 3.8 x 10^^ ergs~^. This 
would require a very special tuning of the astrophysical 
parameters of the system. Yet, even though rare, this 
situation can happen: 4U2206+54 is a low luminosity 
(Lx « 10^5-36 ergs-i) HMXB with a NS accreting from 
the wind of a massive 09.5V star while the X-ray spec- 
trum presents very few (if any) local absorption (Reig 
et al. 2012). Contrary to the case of 4U2206+54, how- 
ever, no pulsations have been found for HDl 19682 or any 
other 7 Cas analog so far. In contrast, within the frame- 
work of the standalone active star, the single absorption 
column, compatible with the ISM value, can be simply 
understood if the system is seen at very high inclina- 
tion angles. It has been recently found that the hard 
X-ray component of the multi-thermal emission of 7 Cas 
(Smith et al. 2012) and HD110432 (Torrejon et al. 2012) 
is very likely produced near the surface of the Be star and 



is absorbed by material with a column density compati- 
ble with the circumstelar disk. If this is also correct for 
HD119682, then the « Njf^^ implies that the sys- 
tem must be observed at high inclination. This would be 
consistent with the single peaked emission lines seen in 
the optical spectrum of this object (SOT). 

R06 and SOT report hints of Fe Ka emission. However, 
no sign of such emission line (nor any other) is present 
here. Thus, the fluorescence emission from neutral Fe 
seems to be variable in the long term. On the other hand 
the 0.5-8.5 keV X-ray luminosity of the source has varied 
within a maximum factor of ~ 2 in a period of seven years 
(R06, SOT, this work). One of the possible sites for the 
Ka reprocessing of the X-rays is the circiimstelar disk. 
Smith et al. (2012) find a direct correlation between the 
absorption column of the hard emission and the EW(Fe 
Ka) (as well as with the reddening) for 7 Cas (Table 5 
of that reference) implying that most, perhaps all, of the 
X-ray emission is produced behind the circumstellar disk 
material. These authors also deduced an inclination of 
i = [43±1]° for 7 Cas. In the case of HD119682, however, 
the absorption column remains always compatible with 
the ISM value despite the changes in the X-ray brightness 
of the source and variable Fe Ka emission because the 
circumstellar disk, confined to the equatorial plane, never 
intercepts the X-ray emission. 

HD119682 would thus appear as a pole-on 7 Cas, an 
important benchmark to explore the phenomenon where 
the plasma, and its variations, could be directly observed 
without any intervening circumstellar material. 
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